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In a recent Letter Shukla and Eliasson reported a 
new attractive force between ions in degenerate quantum 
plasmas. This interesting result has potentially signifi- 
cant implications, as it can lead to formation of ion clus- 
ters in quantum plasmas, even at densities a few times 
below the solid density Q ■ The suggested important con- 
sequence is the modification of the ignition criterion in 
the fast ignition schemes in inertial confinement fusion. 

It is not surprising that this result has already 
prompted a lively discussion in the literature ^3, 4] . One 
of the reasons is that this prediction is a theoretical one, 
with no experimental verification done so far. Moreover, 
the generalized quantum hydrodynamical model used to 
obtain this result ignores some physics that, when taken 
into account, may significantly affect the predicted at- 
tractive force. Namely, this model, as any hydrodynami- 
cal model, completely misses kinetic effects, which in case 
of degenerate electrons are known to lead to the Kohn 
singularity in the dielectric response function, which in 
turn leads to Friedel oscillations in the potential of a 
charge screened by degenerate electrons 0,0]. The depth 
of (at least some of) the potential wells associated with 
the Friedel oscillations could be comparable or even ex- 
ceed the maximum depth of the potential well obtained 
within the generalized hydrodynamical model of Ref. [l[ , 
as the latter is rather shallow. It is thus important to 
compare the potentials produced by both models to es- 
tablish whether and when the predicted attractive part of 
the Shukla-Eliasson's (SE) screened ion potential is im- 
portant relative to the attractive parts of such "Kohn" 
potential predicted by the kinetic model. 

The generalized quantum hydrodynamical model 
(GQHD model) used in Ref. [l[ yields a potential SE 
whose shape depends on the coupling parameter a of the 
electron gas with exchange interaction and correlations 
taken into account [l|: 
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where as — H 2 /e 2 m is the Bohr radius, and ro = 
is the average inter-electron distance in the plasma. The 
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most interesting case is a > 1/4, when the resulting po- 
tential <fi SE has an attractive part. This happens in a 
limited range of as/ro (see Fig. [2]). For a > 1/4, <j) SE 
has the following dependence on the dimensionless dis- 
tance R = r\fiuj pe / vp from the test charge (here ui pe is 



the electron plasma frequency, vf 
the electron Fermi velocity): 
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with A lj2 = (a/a ) 1/4 [y/Aa =p l/±d)^", b = 1/y/Aa - 1, 
and ao is the coupling parameter of an electron gas with 
exchange interactions and correlations neglected (ideal 

1 /3 

electron gas), a a = (9/7r 5 ) (vo/cle)- The attractive 
part of ^„>i/4 is most pronounced when the coupling 
parameter a is at its maximum, a ~ 0.627 at as/ r o ~ 
0.15, corresponding to no ~ 2 x 10 22 cm" 3 , a few times 
below solid densities @. 

The quantum kinetic model (Wigncr model) [B], 0- 
Q accounts for (i) electron degeneracy (via the explicit 
form of Fermi-Dirac distribution of electrons at equi- 
librium), (ii) quantum recoil (explicitly included in the 
Wigncr equation), and (iii) kinetic effects arising from 
the shape of the degenerate electron distribution func- 
tion (e.g;, Kohn singularity in the static response func- 
tion p,Qi|)- Yet this model, in its simplest form, ignores 
electron exchange interaction and correlations. The po- 
tential Kohn of a test charge Q, yielded by this model 
is 
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where K = kvp/VSoj. 
e{0,K) at K 



Due to the Kohn singularity of 



y/3/ao, the potential (j3|) has an oscil- 
latory structure at sufficiently large R, i.e., far from the 
test charge (Friedel oscillations) 

Both potentials © and ^ have attractive parts: 
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FIG. 1: (Color 
(vF/V3uj pe )4>/Q as 



online) The normalized potential 
a function of the normalized dis- 



tance R = r\/ r 3uj pe /vF from the test charge Q, obtained 
from the GQHD model [Eq. (J5|, dashed curve] and from 
the kinetic model [Eq. (|3]), solid curve]. The potentials are 
plotted for as/ro « 0.15, at which the coupling parameter 



([T]) is maximum, 
deepest. 



and the attractive part of 
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is the 



The dependence of rj on / ro is shown in Fig. [2j It is 
seen from Figs [T] and [5] that the depth ratio of the SE 
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FIG. 2: (Color online) The attractive potential well depth 
ratio rj ((HJ (solid green line) and the coupling parameters a 
(dashed blue line) as functions of as/ro. The critical value 
a = 1/4 is indicated with a dotted red line, and the dotted 
gray line marks the unity, for visual reference. 



has a series of potential wells associated with Friedel os- 
cillations (see Fig. [1]). The reasons behind these wells 
are, however, different. The well in ^^^(fl) occurs 
as a result of the interplay between quantum statisti- 
cal pressure and quantum recoil on one hand, and ex- 
change interaction and correlations of electrons on the 
other hand; it vanishes if the latter are neglected. The 
wells in 4> Kohn (R) occur due to the quantum kinetic ef- 
fects that are neglected in the GQHD model. A com- 
parison of the respective attractive well depths of the 
potentials ^ and © is shown in Fig. Q] for the case of 
maximum a (a^/Vo ~ 0.15), when the attractive well of 
<^a>i/4(-R) is the deepest. 

Let us define a parameter characterizing the relative 
importance of the attractive potential wells of <^j,>i/4(-^) 

and <£ Kohn (i?) as v = | A<^ 1/4 (i^ n )/A<£ Kohn (i?- 
wherc A0SE i/4(jR sej and A K ° h "(i? ; 
tivcly, the values of ^> 1/4 (i?) and Kolln (i?) 
bottoms of their respective potential wells at R^ n and 
R^ n , closest to the charge Q. 
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attractive potential well and of the first attractive well 
in the "Kohn" potential is less than unity, for all plasma 
densities including that at which the Shukla-Eliasson's 
attraction is predicted to be the strongest (at the maxi- 
mum of a). 



This suggests that the proper treatment of charge 
shielding should necessarily be a kinetic one, as the 
GQHD treatment does not take into account thc impor- 
tant (and dominant ) part of the attractive force. Another 
reason that warrants a kinetic treatment of the prob- 
lem is that the GQHD model's response is only valid 
for k\p <C 1 [H, where Xf = VF/uj pe is the Thomas- 
Fermi length. Yet in Ref. [l| the test charge potential 
was obtained through integrating its Fourier transform 
4>k = 4:irQ/k 2 e(0,k) over the entire range of k from 
to 00, despite the fact that e(0, k) given by their GQHD 
model is not valid at large k. A fully kinetic treatment 
is free from this limitation, yielding the response that is 
valid for any k. However, a proper way of accounting for 
exchange interaction and correlations within the frame- 
work of kinetic theory is not yet fully developed. Further 
work in this direction is being done and results will be 
reported elsewhere. 



This work was partially supported by the Australian 
Research Council and the Ministry of Education of the 
Russian Federation. The authors thank P.K. Shukla and 
B. Eliasson for useful and inspiring discussions. 



3 



[1] P. K. Shukla and B. Eliasson, Phys. Rev. Lett. 108, 
165007 (2012). 

[2] P. K. Shukla and B. Eliasson, J. Plasma Phys. (under 
consideration) . 

[3] M. Bonitz, E. Pehlke, and T. Schoof, arXiv:1205.4922 

[physics. plasm-ph] (2012). 
[4] P. Shukla, M. Akbari-Moghanjoughi, and B. Eliasson, 

arXiv:1206.3456 [physics.plasm-ph] (2012). 
[5] E. M. Lifshitz and L. P. Pitaevskii, Physical Kinetics 

(Pergamon Press, 1981). 



[6] D. Else, R. Kompaneets, and S. V. Vladimirov, Phys. 
Rev. E 82, 026410 (2011). 

[7] Y. L. Klimontovich, Statistical Physics (Harwood Aca- 
demic Publishers, 1986). 

[8] S. V. Vladimirov and Y. O. Tyshetskiy, Phys. Usp. 54, 
1243 (2011). 

[9] Y. Tyshetskiy, S. V. Vladimirov, and R. Kompaneets, 
Phys. Plasmas 18, 112104 (2011). 
[10] S. V. Vladimirov, Phys. Scr. 49, 625 (1994). 



